Dark radiation candidates after Planck data 
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Recent Cosmic Microwave Background (CMB) results from the Planck satellite, combined with 
previous CMB data and Hubble constant measurements from the Hubble Space Telescope, provide a 
constraint on the effective number of relativistic degrees of freedom 3.62^Q'4g at 95% CL. These new 
measurements of A^^ff provide a unique opportunity to place limits on models containing relativistic 
species at the decoupling epoch. Here we review the bounds or the allowed parameter regions in 
sterile neutrino models, hadronic axion models as well as on extended dark sectors with additional 
light species based on the latest Planck CMB observations. 
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I. INTRODUCTION 

Recent Cosmic Microwave Background (CMB) mea- 
surements from the Planck satellite, combined with other 
cosmological data sets, have provided new constraints on 
the effective number of relativistic degrees of freedom 
iVoff, defined in terms of the energy density of the total 
radiation component as 



Prad 



7/4X4/3 



Pi 



(1) 



where p-^ is the current energy density of the CMB. In 
the standard scenario, the expected value is A'cff = 3.046, 
corresponding to the three active neutrino contribution 
and considering effects related to non-instantaneous neu- 
trino decoupling and QED finite temperature corrections 
to the plasma. Planck data 1], combined with measure- 
ments of the Hubble constant Hq from the Hubble Space 
Telescope (HST) [2| give the constraint TVoff = 3.83±0.54 
at 95% CL. When low multipole polarization measure- 
ments from the Wilkinson Microwave Anisotropy Probe 
(WMAP) 9 year data release Q and high multipole 
CMB data from both the Atacama Cosmology Telescope 
(ACT) 3] and the South Pole Telescope (SPT) d Q are 
added in the analysis, the constraint on iVoff is 3.62^048 
at 95% CL d|. These bounds indicate the presence of an 
extra dark radiation component at the ~ 2.4cr confidence 
level. Different cosmological analyses carried out pre- 
viously to Planck data release including SPT data-only 
have shown a similar evidence [S-Q' see also Refs. |3l.[l0l- 
|28| for constraints on the dark radiation abundances ex- 
ploiting different cosmological scenarios, data sets and/or 
analysis techniques. In addition, the presence of an ex- 
tra dark radiation component will help enormously in the 
agreement on the value of the Hubble constant extracted 
from CMB Planck data and the value of Hq measured by 
the HST team [l|. Even if the discrepancy between the 
CMB and the astrophysical measurements of Hq can be 
alleviated in the context of Hubble bubble models [29] , in 
which we, observers, are living inside a local underden- 
sity, it is mandatory to analyse carefully the constraints 



from Planck data on any physical mechanism which could 
provide AA'off ^ 0.6. 

The simplest scenario to explain the extra dark ra- 
diation AA^cff = Ncs — 3.046 arising from cosmological 
data analyses includes extra sterile neutrino species, since 
there is no fundamental symmetry in nature forcing a def- 
inite number of right-handed (sterile) neutrino species. 
Therefore, sterile neutrinos arc allowed in the Standard 
Model fermion content. However, there are other possi- 
bilities which are as well closely related to minimal ex- 
tensions to the standard model of elementary particles, 
as thermal axions, or extended dark sectors with addi- 
tional relativistic species. New Planck data provide a 
unique opportunity to place limits (or find the favoured 
regions) on the different parameters which describe the 
three models listed above or any other model contain- 
ing new light species, see Ref. [SOf. It is the aim of 
this paper to carefully study these limits. Namely, in 
the case of sterile neutrino models, the constraints on 
Ncs from recent Planck data can set upper bounds on 
the sterile neutrino mixing parameters for sterile neu- 
trino masses m < 0.3 eV, see Ref. [3l| for a recent study. 
We shall focus here on the so-called (3-1-1) neutrino mass 
models 3^|. In the hadronic axion model [3j, ISJ], one 
can explore, as a function of the axion mass rua (being 
TTT-a ^ 0.3 eV) if the axion abundance (parameterized in 
terms of ATVcg) agrees with Planck findings. Models con- 
taining a dark sector with light species that eventually 
decouples from the standard model will also contribute 
to A'cff, as, for instance, asymmetric dark matter models 
(see e.g. Refs. [33, [S^ and references therein), or ex- 
tended weakly-interacting massive particle models (see 
the recent work presented in Ref. [33|)- We will fol- 
low the expressions from Ref. (36[ , in which the authors 
have followed a general approach to describe the dark 
sector structure, including both light and heavy relativis- 
tic degrees of freedom in the dark sector at the time of 
decoupling. While the former correspond to the num- 
ber of degrees of freedom that ultimately constitute the 
dark radiation sector, the latter correspond to relatively 
heavy degrees of freedom that will turn non-relativistic 
and heat the dark radiation fluid. We derive here the 
constraints on the number of light and heavy degrees of 



freedom of the dark sector as a function of its decoupling 
temperature from the standard model sector. 

The paper is organised as follows. Section presents 
the constraints on the sterile neutrino mixing parameters 
in the (3 + 1) neutrino mass models. In Section IlIII we 
briefly review the thermal axion model and illustrate the 
constraints on its mass and its coupling parameter arising 
from Planck measurements on N^ff. Section analyses the 
implications from Planck data on extended dark sectors 
models. Finally, we draw our conclusions in Sec. |Vl 



II. LIGHT STERILE NEUTRINO MODELS 

A number of studies in the literature have been de- 
voted to compute constraints on the light sterile massive 
neutrino thermal abundances [3, [ifl, l38l - |42| . However, 
the extra sterile neutrinos do not necessarily need to fea- 
ture thermal abundances, depending dramatically their 
contribution to the mass-energy density of the universe 
on the flavour mixing processes operating at the decou- 
pling period. Such a study was carried out firstly in 
Ref . [4j| , where the authors computed the constraints on 
the sterile neutrino masses and abundances arising from 
a joint analysis of short baseline oscillation and cosmolog- 
ical data. More recently, the authors of Ref. [31| have 
shown that the constraints on Ncs from recent Planck 
data can set upper bounds on the sterile neutrino mixing 
angles. We benefit here from the approximated expres- 
sions provided in Ref. [43| to explore the constraints 
on the sterile neutrino mixing parameters arising from 
Planck results. The approximate expressions derived in 
Ref. :43,] are valid here, as we are assuming small mixing 
both between the active and heavy sectors and between 
the sterile and light neutrino sectors. In other words, 
if the flavor neutrinos Va, ot = &t^-,t,s (where s refers 
to the fourth sterile neutrino) are related to the massive 
base Vi^ i = 1,2,3,4, through a 4 x 4 unitary matrix 
which U: 



^U„ 



(2) 



we are assuming that \UaA\, \Ujs\ ^ 1, with a = e,fj,,T 
and j = 1,2,3. In this case, sterile neutrinos never reach 
complete thermalization and their abundances are much 
lower than the thermal one. The sterile neutrino con- 
tributes to the energy density of the Universe with [431 ■ 



rish"^ ~ 7x io-^( 



Ar 



'41 



6^2 



9a 
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with 



e, pt, T and Arri^^ is taken here as the squared 



mass of the extra sterile neutrino, assuming mi ~ 0. The 
constants Ca (Ce ^ 0.61 and C^^t- ^ 0.17, respectively) 
are related to the effective potential describing the in- 
teractions of neutrinos with the medium. The constants 
5e — 3.6 and g^ = g^ ~ 2.5 are the coefficients of the 
damping factor. The contribution from the extra sterile 



neutrino to the effective number of relativistic degrees of 
freedom reads: 



A7Veff = 



n,h^ 



7M ^ 
8 V 11^ 



^n^h^ 



(4) 



and therefore, using the Planck measurements of AiVoff 
it is possible to set constraints on the sterile neutrino 
mixing parameters, for a given value of the sterile neu- 
trino mass ?Ti4, provided that nii ^ 0.3 eV. We allow for 
an electron (C/e4 ) and muon {Uf^i) flavor content of the 
sterile neutrino, setting Uri — 0. 

Figure [H left (right) panel, shows the 95% CL con- 
straints on the (|C/e4|j |t^/i4|) plane arising from the 
Planck constraints, A^eff = 3.62^0^8 [N^s = 3.83 ±0.54), 
for two possible values of the sterile neutrino mass, mg = 
0.2 and 0.3 eV. Larger values of the sterile neutrino mass 
will not be relativistic at decoupling and therefore they 
can not be tested exploiting the measured value of A^eff by 
Planck: a full Montecarlo analysis would be needed, anal- 
ysis which will be carried out elsewhere 4J] . Notice that 
the relatively large values of the sterile neutrino mixing 
parameters preferred by short baseline oscillation data in 
(3 -I- 1) models are excluded here for 0.1 < nis ^ 0.3 eV. 
We find Uei < 0.07 and C/^4 < 0.06 at the 95% CL for the 
former range of sterile neutrino masses. For lower sterile 
neutrino masses nis < 0.1 eV, higher mixing parameters 
are allowed, but such a low sterile neutrino mass is highly 
disfavored by oscillation analyses. For instance, the best 
fit point to appearance short baseline data in (3 -f- 1) 
models is found at Am|j^ = 0.15 eV^, being Ue4 = 0.39 
and [7^4 = 0.39 [45]. This region of parameters is, how- 
ever, highly disfavoured by recent Planck measurements. 
Nevertheless one should keep in mind that the analysis 
presented here is in the context of (3 -I- 1) models, which 
have been shown to be inadequate to fit global data sets 
and one should use instead (3 + 2) or (3 + 3) models [45J . 



III. THERMAL AXION MODEL 

Here we first briefly review the origin of axions. Quan- 
tum Chromodynamics (QCD) respects CP symmetry, de- 
spite the existence of a natural, four dimensional, Lorentz 
and gauge invariant operator which violates CP. This CP 
violating-term will induce a non- vanishing neutron dipole 
moment, d„. However, the constraint on the dipole mo- 
ment |d„| < 3 X 10~26 gem |46] requires the CP term 
contribution to be negligible. Why is CP not broken in 
QCD? This is known the so-called strong CP problem. 
The most elegant and promising solution to the strong 
CP problem was provided by Peccei and Quinn |47J |. by 
adding a new global U{1)pq symmetry, which is sponta- 
neously broken at an energy scale /a, generating a new 
spinless particle, the axion. The axion mass is inversely 
proportional to the axion coupling constant fa 



ma 



fa l + R 



0.6 eV 



10^ GeV 

la 



(5) 
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FIG. 1. The left (right) panel illustrates the 95% CL allowed regions in the (|C/e4[, |C^m4|) plane arising from the Planck 
measurements of the effective number of relativistic degrees of freedom A'^off = 3.62lQ'4g (A^cff ~ 3.83 ± 0.54) for two values of 
the sterile neutrino mass. 



where R — 0.553 ± 0.043 is the up-to-down quark 
masses ratio and /^ — 93 MeV is the pion decay con- 
stant. Axions can be produced via thermal or non- 
thermal processes the early universe, providing a pos- 
sible (sub)dominant (hot) dark matter candidate. Here 
we focus on hadronic axion models such as the KSVZ 
model Hi,!!!. 

For axion thermalization purposes, only the axion-pion 
interaction will be relevant. To compute the axion de- 
coupling temperature To we follow the usual freeze out 
condition 



T{Td) = H{Td) . 



The average rate vr 



r = 



TT + a is given by [4S 
1 



1024^5 /2J2Ca.^' 



where 



Ca 



1-R 

3(1 + i?) ' 



(6) 



(7) 



(8) 



is the axion-pion coupling constant 48], and 



/ = n-'T-' 



dxidx2-^-^f{yi)f{y2) 



yiy2 

^ , (s--m2)3(5s-2m2) 
ciuj — 



(9) 



where Ua = {C,^/it'^)T^ is the number density for axions in 
thermal equilibrium, f{y) = l/(e^ — 1) denotes the pion 
distribution function, Xi — \'Pi\/T, yi — Ei/T (z = 1,2), 
s — 2{m'^ +T'^ {yiy2 — xiX2^)) , and we assume a common 
mass for the charged and neutral pious, m^ = 138 MeV. 
We have numerically solved the freeze out equation 
Eq. ([B]), obtaining the axion decoupling temperature 
T]j versus the axion mass rua (or, equivalently, versus 



the axion decay constant fa). From the axion decou- 
pling temperature, we can compute the current axion 
number density, related to the present photon density 



410.5 ±0.5 cm-'^ via 



9*siTo) 
9*s{Td) 



(10) 



where g-^s refers to the number of entropic degrees of 
freedom. At the current temperature, g^^siTo) = 3.91. 
The deviation from the expected value of NcB is 3.046 
due to the presence of a thermal hadronic axion is given 
by 



Pa 



AiV.ff-^ = - - — 



3na 
2 n,y 



(11) 



being Ui, the current neutrino number density. Figure [2] 
illustrates the expected AA^eS as a function of the ther- 
mal axion mass. Axions with masses ma ^ 0.3 eV are 
still relativistic at the decoupling epoch and is precisely 
in this range of values the ones in which CMB N^s mea- 
surements can constrain the hadronic axion model. 

If we assume the A^eff = 3.83 ± 0.54, see the right panel 
of Fig. [21 which corresponds to the value arising from 
the combination of Planck data and HST measurements, 
the thermal axion model is disfavoured at the 2(j CL, 
since axion masses larger than 0.4 eV are excluded by 
cosmology |49l - [5l| while axions with masses nia < 0.4 eV 
do not seem to provide the appropriate amount of dark 
radiation, considering Planck and HST data sets exclu- 
sively. However, when other data sets are also added in 
the analysis, as for instance, WMAP polarization data 
plus high multipole data from both ACT and SPT, Ncs 
turns out to be 3.62lg'^g and therefore the hadronic ax- 
ion model with m^ ^ 0.4 eV is perfectly compatible with 
the value measured of A'off (see the left panel of Fig. ^ . 
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FIG. 2. AAi'cft as a function of the thermal axion mass (in eV). The left (right) panel illustrates the constraint A'eff ~ 3.62lg'4g 
(iVeft= 3.83 ±0.54). 



IV. EXTENDED DARK SECTOR MODELS 

Any model containing a dark sector with relativistic 
degrees of freedom that eventually decouples from the 
standard model sector will contribute to iVcff. An ex- 
ample of these models is the so-called asymmetric dark 
matter scenario, which, in general, contains extra radia- 
tion degrees of freedom produced by the annihilations of 
the thermal dark matter c omp onent. We follow here the 
general approach of Ref. [3g|, in which the dark sector 
containts contains both light {gi) and heavy (g^) rela- 
tivistic degrees of freedom at the temperature of decou- 
pling Ty) from the standard model. For high decoupling 
temperature, Tb > MeV, the contribution to the effective 
number of relativistic degrees of freedom reads [36[ 



A7V, 



eff - 



13.56 [gi + ghf' 
g*s{TD)i of 



(12) 



where gi^siTo) is calculated using the approximated ex- 
pression given in Ref. [52] . If the dark sector decouples 
at lower temperatures (To < MeV), there are two pos- 
sibilities for the couplings of the dark sector with the 
standard model: either the the dark sector couples to 
the electromagnetic plasma or it couples to neutrinos. In 
the second case, which is the one we illustrate here. 



iVe.= (3+^^ 



-gi) 



-)(- 



,3x 



gn+gh 



3 X J + p/j + gt 



gi.i 

_)3 



(13) 
being gH the number of degrees of freedom that become 
non relativistic between Big Bang Nucleosynthesis and 
the dark sector decoupling period. 

As firstly illustrated in Ref. [3^, it is possible to use 
the measured value of N^ff to find the required heavy 
degrees of freedom heating the light dark sector plasma 



gh as a function of the dark sector decoupling temper- 
ature Tjj for a fixed value of ge. Figure [31 left (right) 
panel, illustrates the 2a required ranges for 5^ using 
TVoff = 3.62j:°j° (TVeff = 3.83 ± 0.54), for gn = 0. Notice 
that at decoupling temperatures To > MeV, the stan- 
dard model relativistic degrees of freedom will he heated, 
requiring therefore heating in the dark sector to enhance 
the value of AiVegf. On the other hand, at low decoupling 
temperatures, the number of the required heavy degrees 
of freedom gh decreases as AA^cff does. Indeed, for the 
case of TVcff = 3.62t° j° (TVoff = 3.83±0.54), having extra 
heavy degrees of freedom is highly (mildly) disfavoured. 
This is because at low temperatures, the photon back- 
ground can not get extra heating from standard model 
particles and therefore an extra heating in the dark sector 
will increase dramatically the value of iVeff. 



V. CONCLUSIONS 

Recent Cosmic Microwave Background measurements 
from the Planck satellite, combined with measurements 
of the Hubble constant from the Hubble Space Telescope 
(HST) have provided the constraint N^g = 3.83 ± 0.54 
at 95% CL. If low multipole polarization measurements 
from the Wilkinson Microwave Anisotropy Probe 9 year 
data release and high multipole CMB data from both 
the Atacama Cosmology Telescope and the South Pole 
Telescope are added in the analysis, the constraint on 
TVoff is 3.62lo;4g at 95% CL. These bounds indicate the 
presence of an extra dark radiation component at the 
~ 2(7 confidence level and can be exploited to set limits 
on any model containing extra dark radiation species, 
as sterile neutrino models, hadronic axion scenarios or 
extended dark sector schemes. 

Within the (3-1-1) sterile neutrino scenario, we find that 
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FIG. 3. The left (right) panel shows the 2a required ranges for the number of heavy degrees of freedom heating the dark sector 
Qh using A'eff = 3.62I0 4g (A'off = 3.83 ± 0.54) for several values of gi, the light degrees of freedom of the dark sector. 



the relatively large values of the sterile neutrino mixing 
parameters preferred by short baseline oscillation data in 
(3 + 1) models are excluded here for 0.1 < ttt-s < 0.3 eV. 
For lower sterile neutrino masses rris < 0.1 eV, higher 
mixing parameters are allowed, but such a low sterile 
neutrino mass is highly disfavored by oscillation analy- 
ses. However, other sterile neutrino models, as the (3 + 2) 
or the (3 + 3) scenarios, may provide a much better fit 
to both cosmological measurements and short baseline 
data. In the context of the hadronic axion model, the 
constraint A^cff = 3.83 ±0.54 disfavours the former model 
at the 2cr CL. On the other hand, the axion model studied 
here with ma ^ 0.4 eV is perfectly compatible with cos- 
mological data when lower values of N^s are considered, 
as those obtained when other data sets are analysed to- 
gether with Planck data. Concerning models with a dark 



sector with light species that eventually decouples from 
the standard model, as, for instance, asymmetric dark 
matter models, having extra heavy degrees of freedom in 
the dark sector is highly (mildly) disfavoured for the case 
of TVoff = 3.62to js (^off = 3.83 ± 0.54). Future Planck 
polarization data will help in cornering dark radiation 
models. 
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